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Abstract 

The satellite-borne experiment PAMELA has been used to make a new measurement of the 
cosmic-ray antiproton flux and the antiproton-to-proton flux ratio which extends previously pub- 
lished measurements down to 60 MeV and up to 180 GeV in kinetic energy. During 850 days of 
data acquisition approximately 1500 antiprotons were observed. The measurements are consistent 
with purely secondary production of antiprotons in the galaxy. More precise secondary production 
models are required for a complete interpretation of the results. 

PACS numbers: 96.50.sb, 95.35. +d, 95.55.Vj 



Antiprotons and positrons are a small but not negligible component of the cosmic radia- 
tion. They can be produced in the interactions between cosmic-ray nuclei and the interstellar 
matter. Detailed measurements of the cosmic-ray antiproton energy spectrum therefore pro- 
vide important information concerning the origin and propagation of cosmic-rays. Exotic 
sources of primary antiprotons such as the annihilation of dark matter particles [lH3| and 
the evaporation of primordial black holes [4, |5j can also be probed. The theoretical energy 
spectrum of secondary antiprotons has a distinct peak around 2 GeV and rapidly decreases 
towards lower energies due to the kinematic constraints on the antiproton production. At 
ligher energies the spectrum is slightly steeper than that of the parent protons (e.g. see 
6|), which results in a slight decrease of the antiproton-to-proton flux ratio. 

Since July 2006, PAMELA (a Payload for Antimatter Matter Exploration and Light- 
nuclei Astrophysics) is measuring the antiparticle component of the cosmic radiation. A 
previous PAMELA measurement of the antiproton-to-proton flux ratio between 1.5 and 100 
GeV [7j, was found to follow the expectation from secondary production calculations. How- 
ever, the positron fraction p, |9[ measured in the same energy range showed a clear deviation 
from secondary production models. In order to explain these results both astrophysical ob- 



jects (e.g. pulsars) and dark matter have been proposed as positron sources (e.g. [10]). A 
contribution from pulsars would naturally increase the positron and electron abundances 
without affecting the antiproton component. Other astrophysical models [ll| have been 
proposed to explain the PAMELA positron results but produce an increase in the antipro- 
ton component at very high energies (>100 GeV). A dark matter contribution may require 
pure leptonic annihilation channels, e.g. [12|, or the introduction of a new dark sector of 
forces, e.g. 13|. In [14| it is noted that any signal in the antiproton energy spectrum may be 
hidden due to incomplete modelling of secondary production and cosmic-ray propagation. 
A detailed measurement of the antiproton energy spectrum over a large energy range is 
therefore of great interest. 



The PAMELA experiment 10, Il5] comprises (from top to bottom): a time of flight 
system, a magnetic spectrometer with silicon tracker planes, an anticoincidence system, 
an electromagnetic imaging calorimeter, a shower tail catcher scintillator and a neutron 
detector. These components are housed inside a pressurized container attached to the Rus- 
sian Resurs-DKl satellite, which was launched on June 15 th 2006. The orbit is elliptical 
and semi-polar, with an inclination of 70.0° and an altitude varying between 350 km and 



610 km. 

We report on the cosmic-ray antiproton flux over the widest energy range ever achieved: 
60 MeV to 180 GeV. We also confirm and extend the previously published PAMELA 
antiproton-to-proton flux ratio measurement J7| to the same energy range. Data were 
acquired from July 2006 to December 2008 (850 days), corresponding to > 10 9 triggers. 
Triggered events were selected for analysis if the reconstructed rigidity exceeded the vertical 
geomagnetic cut-off (estimated using the satellite orbital information) by a factor of 1.3. 
Downward-going charge-one particles were selected using the time-of-flight and spectrome- 
ter data. Time-of-flight information was also used to select low velocity (anti)protons while 
electrons were rejected using the electromagnetic calorimeter information, as described in 
[7|] . The remaining electron contamination was estimated to be negligible while contamina- 
tion from locally produced pions was found to be about 10% between 1 and 3 GV/c and 
negligible at lower and higher rigidities [7J, |l6j • 

The highest energy at which antiprotons can be unambiguously measured by PAMELA 
is determined by the contamination of "spillover" protons which are reconstructed with an 
incorrect sign of curvature either due to the finite spectrometer resolution or scattering in the 
spectrometer planes. To reduce this contamination, strict requirements were applied on the 
quality of the tracks reconstructed in the spectrometer. For example, tracks accompanied 
by 5-rsy emission were discarded to avoid poorly reconstructed coordinates on the silicon 
planes of the spectrometer. For each track the maximum detectable rigidity (MDR) was 
evaluated on an event-by-event basis by propagating the estimated coordinate errors and 
taking into account the track topology. The MDR was required to be 6 times larger than 
the measured rigidity. This allowed the antiproton measurement to be extended up to 
180 GV/c with acceptable contamination from spillover protons. The contamination was 
estimated using the GPAMELA detector simulation which is based on the GEANT3 package 
[l7j. The simulation contains an accurate representation of the geometry and performance of 



the PAMELA detectors. For the spectrometer 18j the measured noise of each silicon plane 
and performance variations over the duration of the measurement were accounted for. The 
simulation code was validated by comparing the distributions of several significant variables 
(e.g. coordinate residuals, \ 2 an d the covariance matrix from the track fitting) with those 
obtained from real data. The high-energy region of the deflection distribution was studied 
before applying the MDR selection and agreement within 20% was found between data 



and simulation. This difference was taken as a systematic uncertainty on the spillover 
contamination which was estimated to be ~ 30% for the rigidity interval 100-180 GV/c. 
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he efficiencies were carefully studied using both experimental and simulated data [16, 



20J. The time dependence of the detector performance (and therefore also efficiency) 
was studied using proton samples collected during 2 month long periods. The average global 
selection efficiency was measured to be ~ 30%. The number of (anti)protons rejected by 
the selection criteria due to interactions and energy loss within the detector systems was 
estimated using the simulation. The number of antiprotons lost due to this selection is energy 
dependent and varies from ~ 10% below 1 GeV to ~ 6% above 50 GeV. The antiproton flux 
was obtained by considering the geometrical factor (estimated both analytically and with 
simulations) and the total live time which is provided by an on-board clock that times the 
periods during which the apparatus is waiting for a trigger. 

The energy-binned antiproton fluxes and antiproton-to-proton flux ratios are given in 
Table [H The spectrometer resolution has not been unfolded and a systematic uncertainty 
is included to account for this. Contamination from pions and spillover protons has been 
subtracted from the results. The first and second errors in the table represent the statistical 
and systematic uncertainties, respectively. The total systematic uncertainty was obtained 
quadratically summing the various systematic errors considered: acceptance, contamination, 
efficiency estimation, energy losses, interactions and spectrum unfolding. 

Figure [1] shows the antiproton energy spectrum and Figure [2] shows the antiproton-to- 



proton flux ratio measured by PAMELA along with other recent experimental data 
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25J, 



28] and theoretical calculations assuming pure secondary production of antiprotons during 
the propagation of cosmic rays in the galaxy. The curves were calculated for solar minimum, 
which is appropriate for the PAMELA data taking period, using the force field approximation 



29J, |37|. 



The PAMELA results reproduce the expected peak around 2 GeV in the antiproton 
flux and are in overall agreement with pure secondary calculations. The experimental un- 
certainties are smaller than the spread in the different theoretical curves and, therefore, 
provide important constrains on parameters relevant for secondary production calculations. 



For example, the antiproton flux bands from Donato et al. [26| presented in Figure [T] show 
uncertainties on the propagation parameters (dotted lines) and antiproton production cross- 
sections (dashed lines) and indicate larger uncertainties than those present in the PAMELA 



TABLE I: Summary of antiproton results. Antiproton fluxes (xl0~ 3 particles/(m 2 sr s GeV)) and 
antiproton-to-proton flux ratios (xl0~ 5 ). The upper limits are 90% confidence levels. The first 
and second errors represent the statistical and systematic uncertainties, respectively. 
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FIG. 1: The antiproton energy spectrum at the top of the payload obtained in this work compared 
with contemporary measurements 
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■ 251 ] and theoretical calculations for a pure secondary pro- 
duction of antiprotons during the propagation of cosmic rays in the galaxy The dotted and dashed 
lines indicate the upper and lower limits calculated by Donato et al. 26[ for different diffusion mod- 
els, including uncertainties on propagation parameters and antiproton production cross-sections, 
respectively. The solid line shows the calculation by Ptuskin et al. |27] for the case of a Plain 
Diffusion model. 

measurements. Figure [3] shows the PAMELA antiproton-to-proton flux ratio compared with 
a calculation [14| (dashed line) including both a primary antiproton component from the 
annihilation of 180 GeV wino-like neutralinos and secondary antiprotons. This model, based 
on the non-thermal production of dark matter in the early universe, was proposed to ex- 
plain the high-energy rise in the PAMELA positron fraction 8j. As shown by the dashed 
line in Figure [3j a reasonable choice of GALPROP [31] propagation parameters (dashed- 
dotted line) allows a good description of PAMELA antiproton data with the inclusion of 
the wino-annihilation signal. Given current uncertainties on propagation parameters, this 
primary component cannot be ruled out. It has also been suggested that the PAMELA 
positron data can be explained without invoking a primary component. This is possible if 
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FIG. 2: The antiproton-to-proton flux ratio at the top of the payload obtained in this work 
compared with contemporary measurements 
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281 ] and theoretical calculations for a pure 



secondary production of antiprotons during the propagation of cosmic rays in the galaxy. The 
dashed lines show the upper and lower limits calculated by Simon et al. |6j for the Leaky Box 
Model, while the dotted lines show the limits from Donato et al. 1 301 for a Diffusion Reacceleration 

n 

with Convection model. The solid line shows the calculation by Ptuskin et al. [27( for the case of 
a Plain Diffusion model. 

secondary production takes place in the same region where cosmic rays are being acceler- 



ated [Uj. An increase in the antiproton |32j and secondary nuclei abundances 33] are also 
predicted in this model. The solid line in Figure [3] shows the prediction for the high-energy 
antiproton-to-proton flux ratio. While this theoretical prediction is in good agreement with 
the PAMELA data, in this energy region it does not differ significantly from the expecta- 
tion for standard secondary production models. Comparisons with experimental secondary 
cosmic-ray nuclei data are needed along with higher energy antiproton measurements. New 
data on the boron-to-carbon ratio measured by PAMELA will soon become available, while 
the antiproton spectrum is likely to be probed at higher energies by AMS-02 experiment 34 ] 
which will soon be placed on the International Space Station. 



We have measured the antiproton energy spectrum and the antiproton-to-proton flux 
ratio over the most extended energy range ever achieved and with no atmospheric overbur- 
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FIG. 3: The antiproton-to-proton flux ratio at the top of the payload obtained in this work 
compared with theoretical calculations. The dotted lines show the upper and lower limits calculated 
for a pure secondary production of antiprotons during the propagation of cosmic rays in the galaxy 
by Donato et al. |30l for a Diffusion Reacceleration with Convection model. The dashed line is a 

n 

calculation by Kane et al. [lj] including both a primary antiproton component from annihilation 
of 180 GeV wino-like neutralinos and secondary antiprotons (dashed-dotted line for the secondary 



component). The solid line show the calculation by Blasi and Serpico 



for secondary antiprotons 



including an additional antiproton component produced and accelerated at cosmic-ray sources. 



den. Our results are consistent with pure secondary production of antiprotons during the 
propagation of cosmic rays in the galaxy. We note that the quality of our data surpasses the 
current precision of the theoretical modeling of the cosmic-ray acceleration and propagation 
mechanisms. Improved models are needed to allow the full significance of these experimental 
results to be understood. 
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